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We propose a unified framework that reconciles the stunning success of MOND on galactic scales 
with the triumph of the ACDM model on cosmological scales. This is achieved through the physics 
of superfluidity. Dark matter consists of self-interacting axion-like particles that thermalize and 
condense to form a superfluid in galaxies, with ~mK critical temperature. The superfluid phonons 
mediate a MOND acceleration on baryonic matter. Our framework naturally distinguishes between 
galaxies (where MOND is successful) and galaxy clusters (where MOND is not): dark matter has a 
higher temperature in clusters, and hence is in a mixture of superfluid and normal phase. The rich 
and well-studied physics of superfluidity leads to a number of striking observational signatures. 


The standard A Cold Dark Matter (ACDM) model 
does very well at fitting large scale observables. On galac¬ 
tic scales, however, a number of challenges have emerged. 
Disc galaxies display a tight correlation between total 
baryonic mass and asymptotic velocity, My, ~ v*, known 
as the Baryonic Tully-Fisher Relation (BTFR) |I] [2]. Hy- 
drodynamical simulations can reproduce the BTFR by 
tuning baryonic feedback processes, but their stochastic 
nature naturally results in a much larger scatter [3]. Fur¬ 
thermore, the mass ms and phase-space EHS] distribu¬ 
tions of dwarf satellites in the Local Group are puzzling. 

A radical alternative is Modified Newtonian Dynam¬ 
ics (MOND) [TO], which replaces dark matter (DM) with 
a modification of gravity at low acceleration: a ~ on 
(cin do); a ~ ^aNdo (dN •C do), with best-fit value 
cio — 1.2 x 10 -8 cm/s 2 . This empirical force law has 
been remarkably successful at explaining a wide range of 
galactic phenomena m- In the MOND regime, a test 
particle orbits an isolated source according to v 2 /r = 
^G^M h a 0 /r 2 . This gives a constant asymptotic veloc¬ 
ity, v 2 = y/G^My,ao, which in turn implies the BTFR. 

The empirical success of MOND, however, is limited 
to galaxies. The predicted temperature profile in galaxy 
clusters conflicts with observations jT2|. The Tensor- 
Vector-Scalar (TeVeS) relativistic extension [13] fails to 
reproduce the CMB and matter spectra [Pfl [15] . The 
lensing features of merging clusters PS HI] are problem¬ 
atic [18] . This has motivated various hybrid proposals 
that include both DM and MOND, e.g., [T9ll2.'{| . 

In this Letter, together with a longer companion pa¬ 
per jS], we propose a novel framework that unifies the 
DM and MOND phenomena through the physics of su¬ 
perfluidity. Our central idea is that DM forms a super- 
fluid inside galaxies, with a coherence length of order the 
size of galaxies. The critical temperature is ~ mK, which 
intriguingly is comparable to Bose-Einstein condensation 
(BEC) critical temperatures for cold atom gases. Indeed, 
in many ways our DM behaves like cold dark atoms. 

The superfluid nature of DM dramatically changes its 
macroscopic behavior in galaxies. Instead of evolving 
as independent particles, DM is more aptly described as 
collective excitations. Superfluid phonons, in particu¬ 
lar, mediate a MOND-like force between baryons. Since 


superfluidity only occurs at low enough temperature, 
our framework naturally distinguishes between galaxies 
(where MOND is successful) and galaxy clusters (where 
MOND is not). Due to the larger velocity dispersion in 
clusters, DM has a higher temperature and hence is in a 
mixture of superfluid and normal phases |25H27| . 

The superfluid interpretation makes the non-analytic 
nature of the MOND scalar action more palatable. The 
Unitary Fermi Gas, which has attracted much excitement 
in cold atom physics 128 ] . is also governed by a non- 
analytic kinetic term [ 25 ] ■ Our equation of state P ~ p 3 
suggests that the DM superfluid arises through three- 
body interactions. It would be fascinating to find precise 
cold atom systems with the same equation of state, as 
this would give important insights on the microphysical 
interactions underlying our superfluid. Tantalizingly, this 
might allow laboratory simulations of galactic dynamics. 

The idea of DM BEC has been studied before [50HM] , 
with important differences from our work. In BEC DM 
galactic dynamics are caused by the condensate density 
profile; in our case phonons play a key role in explaining 
the BTFR. Moreover, BEC DM has P ~ p 2 instead of ~ 
p 3 . This implies a much lower sound speed, which puts 
BEC DM in tension with observations [35] , 

DM Condensation: In order for DM particles to con¬ 
dense in galaxies, their de Broglie wavelength AdB ~ 
(mv)- 1 must be larger than the interparticle separation 
t ~ (m/p v ir) 1 ^ 3 - From standard collapse theory, the den¬ 
sity at virialization is p v ; r ~ (l + z v i r ) 3 5.4 x 10 -28 g/cm 3 , 
while the virial velocity is v = 113 M\^ 3 i/1 + z v y T km/s, 
where A /12 = M/10 12 Mq. Thus AdB <; i implies 

m < 2.3 (1 + z vil ) 3 ' 8 A/” 1/4 eV. (1) 

The second condition is that DM thermalizes, with 
temperature set by the virial velocity v. The interac¬ 
tion rate is T ~ Mvp vlI ^-, where J\f ~ ^ is the 

/ vu m 1 m (rnv)' i 

Bose enhancement factor. The rate should be larger than 
the inverse dynamical time tdyn ~ p ’ suc ^ that the 
coherence length will span the halo. This translates into 
a bound on the cross section (with m e v = m/eV): 

ajm k (1 + £vir) -7 ^ 2 to 4 v A/^ 3 52 cm 2 /g . (2) 
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FIG. 2: Numerical solution of Lane-Emden equation 

FIG. 1: Fraction of DM particles m the condensate. 


Later on, we will adopt m = 0.6 eV as a fiducial value. 
For M \2 = 1 and z v - lr = 2, the inequality becomes 
a /to > 0.1 cm 2 /g. The lower end is consistent with cur¬ 
rent constraints [551155] on er /to for self-interacting dark 
matter (SIDM) [39], though these constraints must be 
carefully revisited in the superfluid context. 

The critical temperature, obtained by equipartition 
/cb^c = ]r rnv 2 , is in the mK range: 

T c = 6.5 m~v /3 { 1 + z vil ) 2 mK. (3) 

For 0 < T < T c , the system is a mixture of condensate 
and normal components. The fraction of condensed par¬ 
ticles, 1 — (T/T c ) 3 ^ 2 JJD], is shown in Fig. 0 as a function 
of halo mass assuming z v i r = 0. For m < eV, galaxies 
are almost completely condensed while massive clusters 
have a significant normal component. 

Superfluid Phase: The relevant low-energy degrees of 
freedom of a superfluid are phonons, described by a scalar 
field 9. In the presence of a gravitational potential <F, 
the non-relativistic effective action is £ = P(X), where 
X = 9— m$ — (V0) 2 /2m [29]. We conjecture that DM su¬ 
perfluid phonons are governed by the MOND action [TI] 


This is a polytropic equation of state P ~ p 1+1 / n with 
index n = 1/2. In comparison, BEC DM has P ~ p 2 [52] , 
Including phonons excitations 9 = pt+<f>, the quadratic 
action for 0 is £ quad = {f> 2 - ^(V</>) 2 )- The 

sound speed can be immediately read off: 

c s = y / 2p/m • ( 8 ) 

Using ([T]), we compute the static, spherically- 
symmetric density profile of the DM condensate halo. 
Introducing dimensionless variables p = pfE and r = 

\J 327rGNA 2 m 6 £’ with Po denoting the central density, hy¬ 
drostatic equilibrium implies the Lane-Emden equation, 

= ( 9 ) 

with boundary conditions 5(0) = 1 and S'(0) = 0. The 
numerical solution is shown in Fig. [2] It vanishes at £1 ~ 
2.75, which defines the halo size: R = ^j £i- 

Meanwhile the central density is related to the halo mass 
as 02 Po = ilqfa 15^)1 1 wit h ~ -0.5. Combining 

these results, we can solve for po and R: 


£ = |A(2 mf/ 2 Xj\X\ - a 1 X- i 9p h , (4) 

where A is a mass scale, pb is the baryonic matter density, 
and a is a dimensionless constant. This action should 
only be trusted away from X = 0, as we will see later. 
The matter coupling breaks the shift symmetry at the 
1/Mpi level and is thus technically natural. Remark¬ 

ably 0 is strikingly reminiscent of the Unitary Fermi 
Gas, jCufg(A') ~ X 5 / 2 , which is also non-analytic [29] . 

The phonon action 0 uniquely fixes the properties of 
the condensate through standard thermodynamics. At 
finite chemical potential, 9 = pt , and ignoring <F>, the 
pressure is given by the Lagrangian density: 

P(p) = §A(2 mp) 3/2 . (5) 

This is the grand canonical equation of state P = P(p) 

for the condensate. The number density, n = dP/dp , is 

n = A(2 to) 3/ V /2 . (6) 

Combining these expressions with p = mn, we obtain 

P= ( 7 ) 


po ~ M^ 5 mH /5 A^fv 7 x 10 25 g/cm 3 ; 

R ~ A/i 2 /5 m e - v 6/5 A- 2 4 5 36 kpc, (10) 


where A me y = A/meV. Remarkably, with m ~ eV and 
A ~ meV we obtain DM halos of realistic size. Con¬ 
cretely, as fiducial values we will fix 


m = 0.6 eV; A = 0.2meV. (11) 

This implies R ~ 125 kpc for Mom = 10 12 Af 0 . In gen¬ 
eral, we expect this superfluid core to be surrounded by 
a cloud of DM particles in the normal phase, likely de¬ 
scribed by a Navarro-Frenk-White profile [43] . 

Including Baryons: We now derive the phonon pro¬ 
file 9 = pt + <f>(r) in galaxies, modeling baryons as 
a static, spherically-symmetric source. The equation 
of motion, V • (^/2m\X\ = apb(r)/2M P1 , where 

X = p — m$(r) — (f>' 2 /2m, can be readily integrated: 


y/2m\X\ 


aM h {r) _ 
87rA/pir 2 


( 12 ) 
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There are two branches of solutions, depending on the 
sign of X. We focus on the MOND branch (with X < 0): 

<j>'(r) = \Jrn (^p + \J p 2 + n 2 /m 2 ^ , (13) 

where p = fj, — m$. Indeed, for K/m>/iwe have 

<//(r) ~ y/n(r). (14) 

In this limit the scalar acceleration on baryonic particle 
is a^(r) = aA(f)'/Mp\. Matching to the MOND result 
Q-mond = i/ao^N fixes a in terms of A and a o: 

a 3/2 A = \J a 0 M P1 ~ 0.8 meV . (15) 

For the fiducial value © A = 0.2 meV, we get a ~ 2.5. 

As it stands, the X < 0 solution is unstable. It leads to 
unphysical halos, with growing DM density profiles. The 
instability can be seen by expanding Q in <p = to 

obtain C = sign(X) ip 2 + .... Clearly the kinetic 

term is ghostly for X < 0. (The X > 0 branch is stable 
but does not admit a MOND regime [24|.i 

Since DM in actual galactic halos has T ^ 0, however, 
we expect @ to receive finite-temperature corrections 
in galaxies. At finite sub-critical temperature, the La- 
grangian depends on X and two additional scalars [44) 

B = det 1 / 2 <9 At V ,J d /i V’ J i Y = p + c/) + v- V0, (16) 

where ^d(a:, t), I = 1,2,3, and v are respectively the 
Lagrangian coordinates and velocity of the normal fluid. 
For instance, consider the 2-derivative operator 

AC = M 2 (T)Y 2 ~ M 2 (p + </>) 2 , (17) 


where we have specialized to the rest frame of the normal 
fluid, v = 0. This leaves the static profile (131 unchanged, 
but modifies the quadratic Lagrangian by M 2 (p 2 , restor¬ 
ing stability for sufficiently large M. Specifically this is 
the case for M > eV [24], remarkably of the same order 
as to! By the same token it corrects the condensate pres¬ 
sure by an amount A P = M 2 p 2 , which obliterates the 
unwanted growth in the DM density profile, resulting in 
localized, finite-mass halos. 

Another possibility is the finite-T Lagrangian 


P(X, T) = §A(2to) 3/2 AV|A: - j3(T)Y \. (18) 


where /3 > |, as required for stability [24]. The DM con¬ 
densate pressure is P(/z, T) = |\//3 — lA(2mp) 3 / 2 , which 
is identical to ([5]) modulo the redefinition A —> i J/3 — 1A. 
The halo density profile is therefore identical. 

Including baryons, it is easy to show that (12) becomes 


(j)' 2 + 2m (f~ l) A ^ _ 
yjcj)’ 2 + 2m(/3 — l)/t 0 1 j 


(19) 


At small distances (4>' 2 mp), the solution approxi¬ 
mates the MOND profile <j>' ~ y/R ~ 1/r, and this gives 


the dominant force on a test baryonic particle [H]. At 
large distances (</>' 2 <C to/I), the solution tends to <j>' ~ re, 
but this is subdonrinant to the gravitational acceleration 
due to the DM halo. The transition radius tv delineating 
the MOND regime (r <C r*) from the DM-halo regime 
(r ?v) occurs when re ~ mp. Substituting (121 and 
approximating p by its central value Pq/8A 2 to 5 , we find 


n 


M. 


1/10 ( M b \ 2/5 -8/5.-8/15 

’ ’ 1 1 -* 1 " l eV "meV 


11, b Mqm J 


28 kpc, 


( 20 ) 


where Mn j b = Mb/10 7 * * * 11 M 0 . For a galaxy with Mn i b = 
3 and cosmic DM-baryon ratio M DM /Af b = D DM /O b — 
6, the MOND regime extends to r* ~ 70 kpc. 

Relativistic Completion: It is well-known that a su¬ 
perfluid is described in the weak-coupling regime as a self¬ 
interacting complex scalar field with global t/(l) symme¬ 
try. Consider the relativistic field theory 


C = Hd$| 2 -to 2 |$| 2 


A J (|d$| 2 +m 2 |$| 2 ) 3 
3(A2 + |$|2)6 


( 21 ) 


where A c makes <f> = 0 well-defined. Substituting <f> = 
p e i(e+mt) an( ^ taking the non-relativistic limit, we obtain 


C = -{Xpf + 2mp 2 X 


A 4 ((Vp) 2 -2mp 2 X) 3 
6(A =+p 2 ) B 


( 22 ) 


To leading order in gradients we ignore (Vp) 2 contribu¬ 
tions and integrate out p. In the limit A c this gives 


p 2 ~ A V2m (A 2 ) 1/4 = Ay/2m\X\ . (23) 

Substituting this back into ( [22[ ), we recover the MOND 
action (given by the kinetic part of 0 ) for p^> A c . 

A finite A c implies that MOND is restricted to p > A c , 
i.e., ap > ao- Thus the MOND regime does not 
apply to arbitrarily small accelerations. By choosing A c 
a factor of a few smaller than A, the predicted departure 
from MOND can occur around the acceleration scale of 
the Milky Way dwarf spheroidals, which are well-known 
to pose a challenge for MOND [45H50] . 

Cosmology: Given their mass, our DM particles are 
axion-like. The simplest genesis scenario is through vac¬ 
uum displacement, with DM being generated at a time 
when Hi ~ to ~ eV. This corresponds to a photon tem¬ 
perature of ^TeV, i.e. around the weak scale! The DM 
rapidly reaches thermal equilibrium with itself and be¬ 
comes superfluid, but is decoupled from ordinary matter. 
Naturally DM is much colder cosmologically than in 

collapsed structures. The ratio J , which 

is constant cosmologically, can be evaluated at matter- 
radiation equality using p eq ~ 0.4 eV 4 and v eq = Vi^- ~ 

7 C Y, . The result is Sr — 1CT 28 , which is much colder 

y/mM Pi 1 C 

than the typical range ICC 6 — 10 -2 in galaxies. To ob¬ 

tain an acceptable cosmology, we need A and a to assume 

different values cosmologically. We will denote their cos¬ 
mological values by Aq and op- 
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On a cosmological background, the phonon equation 
of motion that derives from Q is given by 

|((2m) 3 W 2 )=^a 3 Pb . (24) 


Since a 3 pb = const., this can be integrated straightfor¬ 
wardly. The resulting (non-relativistic) energy density 
p = mn = mAo( 2 m) 3 / 2 0 1 / 2 is given by 


p = -t 


aoAo 

, j Tfltp b “h Pdust 5 
M pi 


(25) 


where /?dust ~ a~ 3 . In the matter-dominated era, the 
baryonic contribution ~ p^t redshifts as l/a 3//2 . In order 
for the superfluid to behave as ordinary dust, the second 
term should dominate over the first all the way to the 
present time: pdust ^ ^ mtopb- Substituting the age 
of the universe to = 13.9 x 10 9 yrs and assuming a DM- 
to-baryon ratio of /?dust/pb = 6 , we obtain (with m ~ eV) 

a 0 < 2.4 x 1CT 5 ^ . (26) 


Since 6 ~ 8 m^A n ~ « :i - the phonon velocity increases 
with redshift and inevitably results in a breakdown of the 
non-relativistic approximation, 6 <C m. In other words, 
the superfluid becomes relativistic at sufficiently high 
density. This is consistent with the equation of state (|7| : 
w = ^ = 12 ^ 2 m 6 ■ Demanding that it be non-relativistic 
by matter-radiation equality puts a lower bound on Aq 


A 0 > 0.1 eV. 


(27) 


This is roughly four orders of magnitude larger than the 
fiducial value A = 0.2 eV assumed in galaxies. This 
can be achieved, for instance, if A depends on temper¬ 
ature as A(T) = 1+ka (t° t p /4 i with ka ~ 10 4 . Mean¬ 
while, (26) implies a® 10 “ 4 , which is roughly four 
orders of magnitude smaller than the 0(1) value ob¬ 
tained in galaxies. This can be achieved, for instance, 
if a(T) = ao^ + KaPVTc) 1 / 4 ), with n a ~ 10 4 . Note 
that the scale A' ~ aA appearing in the phonon-baryon 
coupling in Q is nearly temperature-independent. 


Gravitational Lensing: In TeVeS m the complete 
absence of DM requires introducing a time-like vector 
field Ap, as well as a complicated coupling between </>, A M 
and baryons in order to reproduce lensing observations. 
In our case, there is no need to introduce an extra vector, 
as the normal fluid already provides us with a time-like 
vector w/L Moreover, our DM contributes to lensing, so 
we are free to generalize the TeVeS coupling. 

Suppose matter fields couple to the effective metric 


- 9» v - 2 a^()>( 7 g^ v + (1 + 7 )u^u v ), (28) 

with 7=1 corresponding to TeVeS. In the weak-field 
limit, this gives h 00 = — 1 — 2(>I> + aA(j)/Mp{), hij = Sij — 
2(4> + ' yaAcj)/Mpi)Sij , where V 2 >I> = 47tGn (pb + Pdm )• 


Hence the lensing signal arises from a combination of the 
u^Uv disformal coupling and the DM condensate density 
profile. Determining the allowed range of 7 will require a 
detailed study, which is beyond the scope of this paper. 
What is clear is that there should be considerably more 
freedom than in TeVeS. It may even be the case that 
7 = — 1 is allowed, such that the coupling to matter 
would reduce to a simple conformal coupling. 

Observational Implications: We conclude with some 
astrophysical implications of our DM superfluid. 
Vortices: When spun faster than a critical velocity, a su¬ 
perfluid develops vortices. The typical angular velocity 
of halos is well above critical [24] , giving rise to an ar¬ 
ray of DM vortices permeating the disc It will be 
interesting to see whether these vortices can be detected 
through substructure lensing, e.g., with ALMA | 52] . 
Galaxy mergers: A key difference with ACDM is the 
merger rate of galaxies. Applying Landau’s criterion, we 
find two possible outcomes. If the infall velocity zq n f is 
less than the phonon sound speed c s (of order the viral 
velocity [ 22 ), then halos will pass through each other 
with negligible dissipation, resulting in multiple encoun¬ 
ters and a longer merger time. If Ui n f <; c s , however, the 
encounter will excite DM particles out of the condensate, 
resulting in dynamical friction and rapid merger. 

Bullet Cluster: For merging galaxy clusters, the outcome 
also depends on the relative fraction of superfluid vs nor¬ 
mal components in the clusters. For subsonic mergers, 
the superfluid cores should pass through each other with 
negligible friction (consistent with the Bullet Cluster), 
while the normal components should be slowed down by 
self interactions. Remarkably this picture is consistent 
with the lensing map of the Abell 520 “train wreck” [53}- 
[56] . which show lensing peaks coincident with galaxies 
(superfluid components), as well as peaks coincident with 
the X-ray luminosity peaks (normal components). 
Dark-bright solitons: Galaxies in the process of merg¬ 
ing should exhibit interference patterns (so-called dark- 
bright solitons) that have been observed in BECs coun¬ 
terflowing at super-critical velocities m- This can po¬ 
tentially offer an alternative mechanism to generate the 
spectacular shells seen around elliptical galaxies [55] . 
Globular clusters: Globular clusters are well-known to 
contain negligible amount of DM, and as such pose a 
problem for MOND [59]. In our case the presence of a sig¬ 
nificant DM component is necessary for MOND. If what¬ 
ever mechanism responsible for DM removal in ACDM 
is also effective here, our model would predict DM-free 
(and hence MOND-free) globular clusters. 
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